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Abstract 

A strong analogy is found between the evolution of localized disturbances in 
extended chaotic systems and the propagation of fronts separating different 
phases. A condition for the evolution to be controlled by nonlinear mecha- 
nisms is derived on the basis of this relationship. An approximate expression 
for the nonlinear velocity is also determined by extending the concept of Lya- 
punov exponent to growth rate of finite perturbations. 
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In the last years the study of front propagation in spatially extended systems has known 
a renewed interest, due to the relevance of spreading fronts for the emergence of spatial 
structures (patterns) in non-equilibrium systems [![]. In particular, simple reaction-diffusion 
models seem to be appropriate for describing propagation phenomena in different fields, such 
as fluid dynamics, liquid crystals ||, epidemics [Q, chemical reactions, crystal growth ||] 
and biological aggregation ||. Several mathematical models which describe the spreading 
of a disturbance into unstable (or metastable) steady states have been studied in detail, in 
order to uncover the mechanisms underlying the propagation of fronts |6| [9| . 

The main result of these studies can be summerized with reference to the one dimensional 
equation 

u t = u xx + g(u) (1) 

where g(u) G C^fO, 1], g(0) = g{l) = 0. If g > in (0, 1), then u = is an unstable fixed 
point, while u = 1 is a stable one. In this case, any sufficiently localized initial perturbation 
u(x, t = 0) generates a propagating front joining the unstable to the stable state (Fig.l). A 
linear stability analysis shows that the front can have any speed vp larger than a minimal 
value t>L which depends on the behavior of g(u) at u — 0. However, very often velocities 
larger than t>L require special initial conditions to be realized, so that the "physical" speed 
is exactly vp = t>L. In the following we shall call t>L the linear velocity. Whether it is selected 
or not depends on the behavior of g(u) for u > 0. In particular, it has been shown in [TIJ 
that convexity of g(u) is sufficient for vp = t>L- Intuitively, we can say that, if g'(0) > g'(u) 
for all u > 0, the front is "pulled" by the initial growth of u and, otherwise, it is "pushed" 
by the faster growth of finite u |TT| . 

In the present Letter we consider a different problem, namely the propagation of some 
perturbation in a chaotic system (Fig. 2). Thus the front does not separate two different 
phases, since the system is chaotic (and hence unstable) on both sides of the front. More 
precisely, we consider two realizations of a 1-d coupled map lattice (CML) which differ 
only locally in the initial conditions, and we watch the spreading of the relative deviation. 
In spite of the obvious difference with the situation discussed above, we will show that 
there are surprising similarities. In particular, the derivative g'(u = 0) will be replaced 
by the Lyapunov exponent. In order to formulate (heuristically) a condition equivalent to 
g'(0) > g'(u), we will introduce a new indicator of the sensitivity to finite perturbations. We 
shall see that there exists again a minimal velocity t>L, and that the "physical" velocity vp 
can be larger that only if this indicator grows with the perturbation. 

The CML is written as 

= m) (2) 



x? = (l-e)x? + -(xSL 1 +xJ hl ) (3) 

where i and n indicate the discrete space and time variables, and e the diffusive coupling 
parameter. We use periodic boundary conditions on a chain of length L, xf = xf± L . The 
function f{x) is assumed to be a map of some interval into itself. We have chosen this 
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for numerical convenience. We are confident that the basic features can be generalized to 
continuous systems. Disturbances spreading with vp > in CMLs have been observed for 
the first time in [|1^]. Here, instead of considering the map G(x) studied in ||13|| , we shall 
discuss two simpler examples: the "generalized Bernoulli shift" 

f(x) = rx mod 1 (4) 

and the circle map 

f(x) = (x + a) mod 1 . (5) 
The propagation of infinitesimal disturbances is governed by the evolution in tangent space, 



< +1 = /K) 



W+ 2 KVi + M i-i 



(6) 



where /' = df jdx. Instead of considering an initially localized perturbation, we shall first 
refer to a perturbation decaying exponentially for i — > oo, 

< ~ e""* . (7) 

Its temporal growth depends on /i, 

< ~ e AMn "^ . (8) 

The position of the front is defined as the rightmost site where is larger than some 
arbitrarily fixed constant 0(1). This gives for its velocity 

an fi 

For an absolutely unstable system we have A(/i = 0) > , [] so that V{jj) diverges for 
a — > 0. This is intuitively obvious: an almost flat front will appear to move with arbitrarily 



large velocity. On the other hand, it can be shown [|TJj] that, for nearest-neighbour coupling, 
V(fi) — » 1 for fi — » oo. 

We now want to determine the speed wl when the initial perturbation is localized near 
i = and still infinitesimal (the case of finite perturbations will be discussed later). Since 
we expect that any front will have a leading edge where it is infinitesimal and exponentially 
decaying with some exponent /io, we have t>L = V(fio). 



throughout this letter, all Lyapunov exponents are maximal ones, and all perturbations are 
assumed to be typical so that they grow with maximal rate. There exist of course also atypical 
perturbations the growth of which is governed by non- leading Lyapunov exponents 14], but they 
will be neglected. 
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To determine fio and t>L, we need the (maximal) co-mov ing Lyapunov exponent A(v) [15 
For a given v, this gives the local growth rate of a disturbance in a reference frame moving 
with velocity v, u™ ~ e A ^- |ra if z = t>n. The selected front speed is such that a disturbance 
neither grows nor decreases at t>L, i.e. A(vi) = 0. In order to express this in terms of A(/x) 
and fi, we recall that they are related to A(i>) through the Legendre transformation |T6|,|i~Tf 



Hv)=Mfl) -/m ; v= m do) 

djj, djj, 



Therefore, the derivative of V(fi) is directly related to the co-moving exponent, 

dV 1 fdX A\ A(v) 



djj, ji \d/j, ji J ji 
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(11) 



Using A(^l) = 0, we now see that dV jd\i = at a value jio for which v (ji) = Vi, and since 
A(v) is convex (being a Legendre transform), this will be the unique minimum of V(ji). 
Finally, we can write 



Po \ dt* ' (12) 

Thus as long as we can consider a perturbation as infinitesimal, it is the lowest possible 
speed which is selected, which justifies us calling it the "linear velocity" . 

This expression for is identical to that found in Ref. || for the propagation into 
unstable steady states, provided that A(/i) and \i are identified with the complex part of the 
frequency and of the wavevector, respectively. Thus, the relation A = A(/x) plays essentially 
the role of a dispersion relation ]T4]Ji7[ ] . 

Recalling that for closed systems, A(v ) is always a decreasing function (limiting us to 
v > for symmetry reasons) and that A(v = 0) = A(0) ||15|| , we can readily deduce from 
Eq. (0) that t> L is defined if and only if the system is absolutely unstable, i.e. A(0) > 0. 
As can be seen from Fig. 3, V(fi) steadily increases with /i and V(fi — » 0) — » — oo if the 
local dynamics is not chaotic (A(0) < 0). A negative velocity indicates that the perturbation 
regresses instead of propagating: the system is absolutely stable. 

Finally, we consider localized and finite initial perturbations. We call the corresponding 
front velocity vf. Since any front will have an infinitesimal leading edge, we have to expect 
that t>F = V(n*) for some value fi*. It is hard to see how ji* could be smaller than /i , 
whence we just have to distinguish two possibilities: the "linear" (or "pulled") case with 
fi* = /i° and vf = t>L, and the "nonlinear" (or "pushed") case with fi* > fi°, vf > v^. 

In order to see which case is realized in a particular model, we simulate two chaotic 
configurations {x™} and {yf} initially differing in a limited region of the chain (typically 50 
sites in chains of > 1024 sites) and coinciding elsewhere. The front position after n iterations 
is defined as 

R(n) = max{z : \x? - y?\ > 6} . (13) 
where 9 is a preassigned threshold << 1. The front velocity is then defined as 
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Vp 



lim 

n— >oo 



R(n) 



n 



(14) 



We have verified that vp is independent of the amplitude of the initial perturbation Sq and 
of the value of the threshold 9 when they are varied from 10~ 14 to 10 _1 . 

In this way we measured Vp and v l for several CML models and couplings. As expected, 
we found always vp > v^. In most cases, vp = v L (this was found for logistic, cubic and tent 
coupled maps for all tested values of the parameters and of e), but we have also identified a 
class of maps (namely, models (f|), (|5|) and the map G studied in Ref. fl3|l ) where the strict 
inequality vp > is found to hold. The common characteristic of these maps is that f'(x) 
exhibits a narrow peak (or even a 5-singularity). Moreover, in system (|J) [18| and in Ref. 



|T5| a transition between the two above regimes is found upon varying a parameter of the 
map. For map (|5|), such a transition cannot occur since t>L is always zero, the map being 
marginally stable. However, also in this case we can observe a finite vp for a range of a and 
e values. This fact stresses even more that this propagation mechanism is not related to 
local chaoticity, i.e. to sensitive dependence on local and infinitesimal perturbations. The 
unpredictability resulting from the spreading of perturbations does not result here from local 
production of entropy but from entropy transport. 

In order to determine when the nonlinear mechanism is likely to prevail against the 
linear one, we reconsider a heuristic conjecture of van Saarlos [§] for fronts propagating into 
unstable steady and homogeneous states. He observed that vp > only if the local growth 
rate of small but finite perturbations increases with their amplitude. 

In our case the linear local growth rate of perturbations is represented in the limit of 
small coupling e by the Lyapunov exponent of the single map An, which can be defined as 



A 



lim ( lot 

5-*0 \ 



f( x + 6/2)-f(x-S/2) 
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=<log|/'(*)|> 



(15) 



where < . . . > is the average over the invariant measure of the map. If we are interested in 
the evolution of finite disturbances A the average growth-rate will be given by 



1(A) = log 



f(x + A/2)-f(x-A/2) 



A 



=< A(x, A) > 



(16) 



Obviously, limA^o-f(A) = A . Let us first consider map ([|). There, we have 



A{x, A) 



log 



(1-rA) 
A 



, An = log(r) , 
Therefore, the indicator / is given by 



if x G [l/r - A/2, 1/r + A/2] = C(A, r) 
otherwise 



1(A) 



dx i>(x) lot 



[1-rA) 



+ 



/ dx v(x) log(r) 

Jx<£C 



(17) 



where v(x) is the invariant measure. 
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The expression is more compact for the circle map, because there Ao = and the invariant 
measure is flat, so that 

(18) 

This is positive for < A < 1/2, and is an increasing function at small A. Therefore at 
small but finite A we have a positive growth-rate in spite of the stability against infinitesimal 
perturbations. This is always the case when we consider maps like @, @ or G. Conversely, 
for all the other maps we looked at (i.e., logistic and tent maps), we found /(A) < Ao for 
all finite A- values (see Fig. 4). Accordingly, we can conjecture that whenever a nonlinear 
propagation mechanism has been observed, the quantity /(A) is an increasing function at 
small A. If, instead, /(A) < Ao for any A, propagation in the correspondig CML will 
take place with velocity v-p = for any coupling costant e. Nonlinear propagation of 
perturbations can arise only if finite disturbances are, in average, amplified faster than 
infinitesimal ones, i.e. by a factor > exp[Ao] during a single iteration. 

In order to give a quantitative estimate of v F we have to take into account the coupling 
between different sites. We have seen that the linear velocity is the minimum value of 
V(fi) which is obtained from the growth rate A(/i). If vp > t>L, the exponential slope /i* 
of the leading edge is larger than the value /io where V(fi) is minimal. The main effect of 
nonlinearities is to change A(/x) into a function A(/i, A) which coincides with it along the 
leading edge of the front (where A is infinitesimal) but becomes different as A becomes large. 
Our main assumption now is that we have just to replace A(/i) with a suitable average over 
A(/x, A). The average has to be taken over the A range where /(A) > A(0) and which thus 
"pushes" the front. 

The main problem in this assumption is of course that A is a fluctuating quantity. In 
order to make it practically applicable, we have to resort to a mean field approximation. 

By assuming that the perturbation decays exponentially as 

A™ = e"^$™ , (19) 

from Eqs.(2) and (3), we find that it evolves in time according to 

A? +1 = \m + IK) - m - ~ A?) I = A? e^"< A ") (20) 

where 

A" = e-* ((1 - e)S? + £ -m-i + $r + i e^)) . (21) 

We now introduce a mean field approximation by assuming that $™ is independent of i, and 
A(x, A) equal to its average over x. This allows us to rewrite Eq. ( p0|) as 

$ n+i = $ n [(1 _ £)+£ cosh^)] e /(A) (22) 

Performing an average over the range D of A where /(A) > 0, we obtain an effective 
Lyapunov exponent 



/(A) = A log 



:i-a) 

A 



6 



A c (//) = log [(1 -e) + ecosh(/i)] + -L [ dA 1(A) , \D\ = [ dA (23) 

and, in analogy with the linear case, 

VM = ^ • (24) 

Just like V(fi), V c (fi) is a convex function with a unique minimum. It is thus natural to 
assume that the selected velocity for the front will be given by the minimum of (p4[) 



v F = min M V c (/x) . (25) 



The value \l c where V c (fi) is minimal would be equal to fi* if Eq. (|T9| ) would hold with the 
same \x in the leading edge and in the pushing region. This, however, need not be the case 
and we indeed found /x c < /x* in general. 

In Fig. 5, the numerical results are reported toghether with the predictions obtained 
from Eq. fl25|) for the circle map with two different values of a. The agreement between 
simulation and theoretical results is reasonably good for large coupling e. However, it can 
be seen that the front propagates only for e larger than a certain threshold e c (a). Equation 
(p5|) does not predict such a transition which can be attributed to the particular structure of 
the invariant measure for model (|5]) for e < e c (a). The invariant measure becomes extremely 
irregular below threshold and this does not allow any more a "synchronization" of the motion 
of the disturbances, as necessary to observe a front propagation. Obviously, this cannot be 
recovered from a mean field analysis. An analogous comparison for map @) with e — 1/3 is 
reported in Fig. 6. The overall behaviour of the velocity provided by Eq. ( f25|) is in agreement 
with that of the measured vf. More precisely, the theoretical predictions are larger than t>L 
for any value of the parameter r. 

In conclusion, we have demonstrated that the propagation of perturbations in chaotic 
systems is very similar to the propagation of fronts between steady states. This includes 
the possibility of "nonlinear" selection of velocity. We have verified that an extremely crude 
estimate of the influence of nonlinearities on the velocity gives surprisingly good agreement 
with simulations of several coupled map lattices. 
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FIGURES 



Fig.l: A typical front connecting an unstable with a stable region. 

Fig. 2: A typical chaotic state (x), a perturbed state (y), and their difference. The front 
separates the perturbed from the not yet perturbed region. 

Fig. 3: Velocities V(/i) versus fi for the coupled piecewise linear maps (§]) with e = 1/3. 
The solid curve refers to the absolutely unstable situation (r > 1), the dashed line to the 
marginally stable case (r = 1) and the dash-dotted one to the absolutely stable case (r < 1). 

Fig. 4: Nonlinearity indicator 1(A) for the single maps: logistic map at the crisis (solid line); 
tent map (dashed line); circle map @ with a = [1 — (y/E — l)/2] (dotted line); generalized 
Bernoulli shift (H) with r = 1.10 (dash-dotted line). 

Fig. 5: Front velocities for circle coupled maps as a function of the coupling parameter e: 
measured velocity vp (crosses) and theoretical prediction vt (circles). Figure (a) refers to a 
map parameter a = [1 — (y/E — l)/2] and (b) to a = [1 — (v5 — l)/2]/8. 

Fig. 6: As in Fig. 5 for coupled piecewise linear map with r > 1 (e = 1/3). In this case the 
linear velocity t>L is reported too (solid line), since it is positive. 
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